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INTRODUCTION 

One  of  the  primary  uses  for  wind  tunnels  is  to  be  able  to  place 
scale  models  of  vehicles  into  the  tunnel  in  order  to  measure  the 
aerodynamic  loads  on  the  vehicle.  This  concept  of  testing  small- 
scale  vehicles  has  been  done  ever  since  the  development  of  the  first 
airplane.  Early  systems  utilized  strain-gage  beam  balances  mounted 
outside  of  both  the  model  and  the  tunnel.  Later  improvements  in 
instrumentation  made  it  possible  to  place  a cantilever  balance 
instrumented  with  strain  gages  inside  the  model. 

The  primary  vehicles  tested  at  the  Naval  Surface  Weapons  Center, 
White  Oak  Laboratory  (formerly  the  Naval  Ordnance  Laboratory) , have 
been  weapon  configurations  that  had  relatively  small  out-of-plane 
loads  in  comparison  to  those  in  the  plane  of  the  angle  of  attack. 
Consequently,  it  was  possible  to  utilize  a system  whereby  the 
interactions  between  the  in-plane  and  out-of-plane  loads  were  nulled 
by  electrically  shunting  appropriate  arms  of  the  Wheatstone  bridge 
circuits  (Ref.  (i)).  This  is  accomplished  by  loading  the  balance  in 
one  plane  and  then  adjusting  the  resistances  in  the  appropriate  legs 
of  the  bridges  in  the  gages  of  the  orthogonal  plane  so  that  there  is 
no  indication  of  a load  in  this  orthogonal  plane.  This  technique 
effectively  eliminated  all  first-order  interactions  and  proved  to  be 
quite  accurate  and  acceptable. 

In  recent  years,  with  the  introduction  of  re-entry  vehicles  and 
quite  small  tolerances  on  the  measurement  of  in-plane  and  out-of-plane 
loads  it  became  necessary  to  improve  the  accuracy  of  the  load 
measurements.  The  scheme  chosen  to  do  this  was  one  that  has  been 
utilized  for  some  time  by  other  tunnel  complexes  where  a wider 
variety  of  vehicle  types  have  been  tested  (Refs.  (2)  — (5) ) . The 


^Shantz,  L,  Gilbert,  B._  D..,  White,  C.E.,"NOL  Wind-Tunnel  Internal 
Strain-Gage  Balance  System,"  NSWC/WOL,  NAVOP.D  Rpt.  2972,  Sep  1953 

(2) 

Hansen,  Raymond  M. , "Evaluation  and  Calibration  of  Wire-Strain- 
Gage  Wind  Tunnel  Balances  Under  Load,"  AGARD  Rpt.  13,  Feb  1956 

^"Calibration  and  Evaluation  of  Multicomponent  Strain-Gage 

Balances,"  Prepared  for  presentation  to  the  NASA  Interlaboratory 
Force  Measurements  Group  meeting  held  at  JPL  16-17  Apr  1974 

(4) 

Smith,  David  L. , "An  Efficient  Algorithm  Using  Matrix  Methods  to 
Solve  Wind-Tunnel  Force-Balance  Equations,"  Langley  Research 
Center,  NASA  TN  D-6860,  Aug  1972 

^Hurlburt,  R.  T.  and  Berg,  D.  E.,  "A  Non-Iterative  Method  of 
Reducing  Wind  Tunnel  Moment  Balance  Data,"  Sandia  Laboratories, 
Albuquerque,  N.  M. , SC-RR-720665,  Mar  1973 
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basic  scheme  reported  on  in  Reference  (4)  has  been  modified  and 
adapted  to  our  needs.  Basically  this  involves  the  use  of  second- 
order  equations  to  fit  for  the  interactions  rather  than  the  electrical 
nulling  of  only  the  first-order  interactions. 

STRAIN-GAGE  BALANCE 

Forces  and  moments  can  be  measured  by  equating  a change  in 
strain  at  a point  on  a cantilever  beam  to  the  force  and/or  moment 
that  caused  *-he  strain.  The  ideal  situation  is  to  construct  a 
balance  so  that  the  tension  or  compression  in  an  individual  section 
is  sensitive  to  only  one  type  of  loading.  It  is  then  possible  to 
locate  a strain  gage  on  the  surface  of  the  balance  and  calibrate  the 
change  in  resistance  of  the  gage  with  respect  to  the  load  that 
caused  it. 

Rather  than  use  one  gage  it  is  better  to  utilize  four  gages 
arranged  in  a Wheatstone  bridge  circuit  as  shown  in  Figure  1. 


A 


This  design  makes  it  possible  to  measure  the  absolute  change  in 
voltage,  makes  the  output  proportional  to  four  times  the  output 
of  one  gage,  heJ.ps  to  minimize  temperature  effects,  and  helps 
to  isolate  the  load.  If  the  four  gages  in  the  bridge  are 
located  on  a balance  section  so  that  the  load  causes  the 
surface  under  gages  1 and  3 to  be  in  tension  (increase  in  gage 
resistance)  and  the  surface  under  gages  2 and  4 to  be  in  compression 
(decrease  in  gage  resistance)  signal  lead  A would  be  at  a lower 


4 


NSWC/WOL/TR  76-26 


potential  and  signal  B at  a higher  potential.  By  means  of  a cali- 
bration this  change  in  potential  can  be  equated  to  the  load  that 
caused  it. 

In  order  to  measure  a force  and  the  moment  it  causes  about  some 
point  on  the  balance,  it  is  necessary  to  measure  the  strain  in  the 
balance  at  two  locations.  For  instance,  in  Figure  2,  if  a force  is 
applied  at  any  point  along  the  beany  it  generates  a moment  which 
varies  along  the  beam. 


FIG.  2 CANTILEVER  BEAM  BALANCE 


The  moments  at  gage  sections  one  and  two  are  equal  to. 


M, 


Ff, 


M, 


FI. 


These  moments  can  be  measured  with  Wheatstone  bridges  located  at 
sections  one  and  two.  Assuming  that  these  are  measured,  the 
moments  can  be  added  to  get 


M1  + m2  = F(£1  + l2) 


(1) 


or  they  can  be  subtracted  to  get 


M1  - M2  = F(£x  - £2) 


(2) 
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Since, 


4^  - l 2 * A*  * known  distance  between  the  gage  sections, 


the  equations  can  be  written  as 


M1  + M2  = F ( AX  4 2 l2) 
Mx  - M2  = FAX 


These  can  then  be  solved  for  the  two  unknowns,  F and 


F 


l 


2 


M,  4-  M, 
x 2 

2F 


AX 

2 


(3) 

(4) 


(5) 

(6) 


Knowing  these,  the  force  and  its  point  of  application,  the  moment 
about  any  point  on  the  balance  can  be  calculated. 

Utilizing  the  above  concepts,  it  is  possible  to  instrument  the 
balance  in  both  the  pitch  and  yaw  planes.  Since  both  planes  are 
instrumented  the  same,  only  the  pitch  plane  will  be  shown.  There  are 
two  different  wiring  circuits  in  use  at  the  White  Oak  Laboratory. 

In  the  first,  the  two-moment  method,  a forward  balance  section  and 
an  aft  section  are  wired  independently  as  shown  in  Figure  3.  By 
treating  the  fore  and  aft  sections  separately  it  is  possible  to 
minimize  the  effects  of  a temperature  gradient  along  the  balance. 

In  this  case  it  is  necessary  to  add  and  subtract  the  moments 
measured  at  each  section  for  use  in  Equations  (5)  and  (6) . This  is 
done  in  the  data-reduction  program. 

In  the  second  wiring  circuit,  the  force-moment  method,  the 
moments  are  added  and  subtracted  electrically  on  the  balance  as  is 
shown  in  Figure  4.  This  has  the  advantage  in  that  the  output  signals 
are  directly  proportional  to  a force  and  a moment  but  has  the  dis- 
advantage in  that  a temperature  gradient  along  the  balance  causes  a 
change  in  the  output  that  is  hard  to  adjust  or  account  for  in  a 
rational  manner. 

Similar  techniques  are  used  to  measure  the  roll  moment  and  axial 
force,  although  the  mechanical  design  of  the  gage  mounting  sections 
becomes  more  intricate.  As  shown  in  Figure  5,  the  roll  section 
usually  consists  of  a cruciform  section,  and  the  axial  sectic  > of  a 
series  of  thin  webs. 

As  mentioned  earlier,  the  objective  of  the  balance  designer  is 
to  design  gage  sections  that  will  respond  to  only  one  type  of  load 
with  no  interactions  between  loads.  If  this  could  be  accomplished, 
a large  part  of  the  rest  of  this  report  would  not  be  necessary;  but, 
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FORWARD  FITCH 
BRIDGE  SIGNAL 


AFT  FITCH 
BRIDGE  SIGNAL 


FORWARD  PITCH 
BRIDGE  SIGNAL 


AFT  PITCH 
BRIDGE  SIGNAL 


FIG.  3 TWO  MOMENT  FITCH  GAGE  SECTIONS 


7 


FIG.  4 


1 


2 


PITCH 

MOMENT 

SIGNAL 


0 


FORCE 


0> 


FIG.  5 ROLL 


AXIAL 

SIGNAL 


AXIAL  SIGNAL 


iD  AXIAL  SECTIONS 
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since  there  are  always  some  interactions  (although  they  are  usually 
very  small)  it  is  necessary  to  account  for  them  in  the  resolution  of 
the  loads  being  measured. 

GENERAL  LOAD  RESOLUTION  TECHNIQUE 

The  general  scheme  in  using  a strain-gage  balance  to  measure 
aerodynamic  loads  consists  of  first  loading  the  balance  with  known 
loads  and  then  relating  its  response  to  these  known  loads  through  a 
set  of  calibration  constants.  These  constants  can  then  be  used  to 
determine  an  unknown  load  applied  to  the  balance  by  going  through 
the  inverse  process.  Since  there  are  six  degrees  of  freedom  in  the 
model's  system,  it  is  necessary  to  use  a six-component  balance  to 
determine  all  of  them.  These  six  loads  are: 


F.,  = normal  force 
N 

Fy  = yaw  force 

My  = pitching  moment 

Mz  = yawing  moment 

Mx  = rolling  moment 

F.  = axial  force 
A 

The  standard  model  axes  and  loads  are  shown  in  Figure  6. 


If  the  assumption  is  made  that  the  variation  of  the  balance  out- 
put with  load  is  only  slightly  nonlinear,  then  only  second-order 
load  effects  need  to  be  calibrated  (this  assumption  is  checked  during 
each  calibration) . All  possible  second-order  load  combinations  are 
as  follows: 


F 

N 

fn 

fy 

fnfy 

F2 

Y 

“y 

fnmy 

fy”y 

«? 

Mz 

fnmz 

FYMZ 

Vz  MZ 

Mx 

F M 
NX 

fymx 

¥x  KzMx 

fa 

F F 
N A 

fyfa 

Va  mzpa 

For  bookkeeping  purposes,  the  27  load  combinations  have  been  assigned 
the  following  subscripts: 
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Load 

Subscript 

F 

N 

1 

fy 

2 

“y 

3 

Mz 

4 

"x 

5 

FA 

6 

p2 

•7 

fn 

1 

Vy 

8 

Vy 

9 

Vz 

10 

F M 
NX 

11 

fnfa 

12 

F^ 
r Y 

13 

fymy 

14 

Load 

Subscript 

fymz 

15 

fymx 

16 

fyfa 

17 

18 

Vz 

19 

Vx 

20 

Va 

21 

4 

22 

MzMx 

23 

M F 
Z A 

24 

4 

25 

Va 

26 

4 

27 

According  to  the  assumption  that  the  highest  order  of  non- 
linearity in  the  balance  output  is  second  order,  the  output  of  the 
ith  channel  can  be  represented  as 
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6.  = k.  .F„  + k.  _F  + 

1 1,1  n 1,2  y 


+ ki,14FYMY  + 


+ k-  0-,F* 

1,27  A 


(7) 


where  k.  . (i  = 1,2,... 6),  (j  = 1,2,... 27)  are  balance  coefficients 
1 • J 

to  be  determined  in  the  calibration  of  the  balance.  Equation  (7) 

can  thus  be  normalized  by  dividing  both  sides  by  k.  . and  then 

1,1 

expressing  it  as  follows. 


where 


= 1 when 


K.  .F„  + K.  ~F„  + 
1,1  N 1,2  Y 


+ K.  0_F" 
1,27  A 


normalized  intei  ":tion  coefficient  when  i ? j 

i = j 


(8) 


til 

M . = 1/k.  . = i component  sensitivity 
11,1 

(L^)  = = uncorrected  load  on  the  i*"*1  component 

For  the  case  where  i = 1,  Equation  (8)  becomes 

<VU  ' FN  + K1,2FY  + •••  + K1,27FA  ' (9) 

Equation  (9)  can  then  be  arranged  as 

FN  ■ <Fn'u  - <K1,2PY  " •••  + K1,27FA>  ■ (10) 

Equation  (8)  consists  of  six  equations  similar  to  Equation  (10)  which 
present  the  true  loads  as  a function  of  the  uncorrected,  indicated 
loads  and  all  of  the  first-  and  second-order  interactions.  These 
equations  are: 


FN  = (V.  " (K1,2FY  + + Kl,27FA) 


u 


F = 
Y 


(Fy)  - (K2,ifn  + K2,3MY  + •••  + K2,27FA) 


u 


My  (MY}..  ” (K3,lFN  + •••  + K3,27FA) 


u 


MZ  = <MZ>  - (K4,1FN  + •••  + K4,27FA> 
u 


“ (Mx\.  ■ (K5,1FN  + •••  + K5,27FA} 


u 


(11) 

(12) 

(13) 

(14) 

(15) 
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' t?  \ _ 

'*A' 

u 


(K6,1FN  + 


+ K6,27F1> 


(16) 


The  remaining  problem  is  one  of  how  to  utilize  these  equations 
in  order  to  obtain  the  true  loads  acting  on  a model  during  a test. 
Since  all  of  the  equations  are  coupled  through  the  first-  and  second- 
order  Interactions,  an  iterative  approach  has  been  utilized  to 
solve  them.  It  is  not  absolutely  necessary  to  use  an  iterative 
scheme  (see  Ref.  (5))  but  in  order  to  handle  a wide  range  of  balances 
for  which  the  relative  magnitude  of  the  nonlinear  terms  varies,  it 
is  the  most  accurate  and  quickest  way  of  solving  them. 


Matrix  relationships  were  used  in  order  to  expedite  the  handling 
of  Equations  (11) -(16)  by  making  the  following  definitions. 


Output 

channel 

matrix 


® * 


[6.]  i = 1,2, ...6 


Load 

matrix 


Y 

MZ 


[L^]  i — 1,2, ...6 


Second  order 
force  and  moment 
product  matrix 


P 


N 

fnfy 

fnmy 


- LiL: 


i s 1,2,.. .6 
j = i,i+l, . . .6 
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Sensitivity 

matrix 


M 


"M1  1 0 0 0 

0 u,  ~ 0 0 

0 o2'2  Pj  , 0 

0 0 0 " y4 , 4 

0 0 0 0 


0 

0 

0 

0 

,J5,5 


0 

0 

0 

0 

0 


6,6 


l 

Kl,2 

Kl,3 

Kl,4 

K1.5 

Kl,6 

First-order 

interaction 

Cl  = 

K2 , 1 
K3,l 

1 

K3,2 

K2,3 

1 

*2.4 

K3,4 

K2,5 

K3,5 

K2,6 

K3.6 

coeff.  matrix 

X 

K4,l 

K4,2 

K4,3 

1 

K4 , 5 

*4.6 

K5 , 1 

K5,2 

K5,3 

K5,4 

1 

*5.6 

Ki 

K6,2 

K6,3 

K6,4 

K6.5 

1 

[*1.7 

K1.8 

Kl,9 

K1.10 

*1,27 

K2 , 7 

K2,8 

K2 , 9 

K2 , 10 

K2,27 

Second-order 

interaction 

C7  = 

K3 , 7 

K3,8 

K4,8 

K3 , 9 
K4,9 

K3,10 

K3,37 

K4,27 

coeff.  matrix 

L 

K4,7 

4,10 

K5,7 

K5,8 

K5 , 9 

K5,10 

K5 , 27 

_^6 , 7 

K6 , 8 

K6,9 

K6 , 10 

K6,27 

Equations  (11)- (16)  can  then  be  written  as. 


L = C.L  + C-P 

(17) 

u 1 2 

or 

V © = C^L  + C2P 

(18) 

Now  if  C2P  is 

subtracted  from  each  side  and  then  each  side  is 

premultiplied 

by  the  inverse  of  C^,  Equation  (18)  becomes, 

C^U®-  C^1C2P  = C'1C1L  = L 

(19) 

or 

L = c"1^®-  C^1C2P 

(20) 

Letting  - u(h)  , the  uncorrected  loads  which  can  always  be 

calculated  directly  from  the  balance  output  channels,  and 
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then 


M = C"1C2 


L = C'1^  - MP 


(21) 


The  values  for  C~*  and  M are  known  from  the  balance  calibration 

(see  next  section) . The  values  for  Lu  can  be  calculated  from  the 

balance  output  channel  readings  and,  for  the  first  approximation, 

P can  be  calculated  using  the  uncorrected  loads,  Lu;  therefore,  for 

the  first  iteration 


L.  = c7'L  - MP 
1 1 u 


where 


P = (L±)  (L.) 

u -1  u 


l 1 , 2 , . . . 6 

j = i . . . 6 


For  the  second  iteration  use 


and 


P E P = (L.)  (L.) 

1 1 1 3 1 


L2  - C'l\  - MP1 


The  n*"*1  iteration  would  then  be 


L = C7LL  - MP  . 
n 1 u n-1 


(22) 

(23) 

(24) 

(25) 

(26) 


These  iterations  should  be  continued  until  the  differences  in  the 
interaction  loads  are  equal  to  the  resolution  of  the  system,  or 

|MP  , - MP  J < | resolution | . (27) 

n-i  n-4 


Generally  the  digital  resolution  of  the  recording  system  is  expressed 
in  counts,  with  one  count  being  the  smallest  unit  of  resolution,  and 
the  sensitivity  of  the  output  channel,  vk,  is  expressed  as  load/ 

count.  In  that  case  the  resolution  is  equal  to  the  sensitivity  ik. 

The  following  limit  has  proven  to  be  adequate  to  obtain  accurate 
results : 


MP 

n 


2 


(28) 


Usually  only  two  or  three  iterations  are  needed  to  obtain  this  limit. 

Some  of  the  computer  routines  used  to  perform  the  work  reported 
on  in  this  report  are  enclosed  in  Appendices  A,  B and  C.  Appendix  C 
contains  the  data-reduction  subroutines  which  are  used  to 
accomplish  the  above  technique. 
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BALANCE  CALIBRATION 


The  purpose  of  conducting  a balance  calibration  is  to  measure  the 
balance  response  to  known  loads  in  such  a manner  that  the  inverse 
procedure  can  be  utilized  to  measure  unknown  loads  as  a function  of 
the  balance  response  to  these  loads.  As  mentioned  earlier,  a second- 
order,  combined  .oad  calibration  has  always  proved  sufficient.  This 
type  of  calibration  is  performed  by  first  determining  the  balance 
response  to  pure  primary  loads  and  then,  to  combinations  of  these 
loads.  This  is  accomplished  by  attaching  a lightweight,  rigid 
calibrating  bar  to  the  balance  and  then  applying  leads  to  the  bar. 

The  purpose  of  the  bar  is  to  provide  a way  of  transferring  loads  to 
the  balance  in  the  same  manner  as  they  are  transferred  f^om  the 
model  to  the  balance.  Primary  loads  are  applied  one  at  a time  and 
the  balance  response  on  all  channels  is  recorded.  As  an  example,  if 
pure  normal  forces  are  applied  to  the  balance  center  by  hanging 
weighto  at  the  center,  all  of  the  loads  other  than  normal  force  drop 
out  of  Equations  (11)- (16)  and  the  remaining  terms  can  be  utilized 
to  solve  for  the  primary  sensitivity  of  the  balance  to  normal  force. 


For  n weights  there  are  n 


sets  of  the  following  equations: 


FN  = M191  “ {K1,7FN)  = wei9ht  {29) 
FY  * “292  - (K2,1FN  + K2,7FN>  ■ 0 ,30) 
«Y  ' 113e3  - (K3,1FN  + K3.7Fn»  - 0 <31) 
MZ  * u464  ' (K4,1FN  + K4,7FN»  * 0 (33) 
MX  - “5S5  ’ (K5,1FN  + K5,7Fn'  - 0 (33) 


Equation  (29)  can  be  solved  by  fitting  the  values  of  as  a function 

of  applied  weight  in  order  to  determine  u ^ and  The  results 

expressed  on  the  other  channels  as  a function  of  normal  force  are 

recorded  and  saved  for  later  reduction.  Each  of  the  other  primary 

loads  can  be  added  individually  and  each  primary  sensitivity 

obtained  in  the  same  fashion.  Once  the  primary  sensitivities  are 

known,  it  is  possible  to  go  back  and  reduce  the  data  recorded  on 

the  other  balance  channels.  For  example.  Equations  (11) -(16)  can 

now  be  fitted  in  order  to  obtain  the  remaining  unknown  coefficients, 

K.  . and  K.  _ (i  = 2 -*■  6)  . After  this  has  been  accomplished,  all  of 
i , / 

the  primary  sensitivities  and  all  of  the  first-order  interaction 
coefficients  are  known. 
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In  the  determination  of  the  second-order  interaction  coefficients 
it  is  necessary  to  apply  combined  loads  to  the  calibrating  bar.  For 
example,  if  both  a normal  force  and  a yaw  force  are  applied 
simultaneously,  substitution  of  the  loads  into  Equations  (11) -(16) 
would  give  the  following  equations: 


F = 
N 


Vi 


F = 
Y 

0 = 


p292 

p393 


0 = 


u696 


(K1,2PY  + K1,7FN  + K1,8FNFY  + K1.13F^> 

<K2,1FN  + K2,7FN  + K2,8FNFY  + K2,13FY) 

(,t3,lFN  + K3,2FY  + K3,7FN  + K3,8FNFY  + K3.13FY> 


(lt6,lFN  + K6,2FY  + K6,7FN  + K6,8FNFY  + K6,13FY) 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


Since  the  terms  K.  Q(i  = 1 - 6)  are  the  only  unknown  terms,  they  can 

1,0 

be  easily  solved  for  as  follows: 


Ki,8 


Vi  - Ki.lFN  - Ki,2FY  - 


F F 
NY 


Ki , 7FN  ~ Ki,13FY 


(41) 


for  i = 1 - 6. 

In  principle,  only  one  combined  load  has  to  be  loaded  in  order  to 
solve  for  K.  Q,  but  in  order  to  minimize  errors  and  to  help  check 

the  assumption  that  second  order  is  the  highest  order  of  nonlinearity, 

several  are  loaded  and  the  series  are  fit  to  obtain  K.  Q(i  = 1 - 6). 

1,0 

The  remaining  second-order  coefficients  are  obtained  in  a similar 
fashion . 

In  order  to  expedite  the  handling  of  the  calibration  data,  the 
following  definitions  were  made. 

60  . = balance  output 

*,i,n  * 

L.  . = loads  applied 

K.  g = balance  coefficients 

X f X, 

where 

l = group  number  (1-27)  the  same  as  used  to  define  the 
27-load  configurations  on  page  12. 

i = output  chani  el  number  (1-6) 

n = sequential  numerical  order  of  the  weights  (1-15) 
t = type  of  load 
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The  correspondence  between  these  values  is  as  follows: 


load 

F 

rN 

fy 

"y 

Mz 

*x 

F. 


Type, 


1 

2 

3 

4 

5 

6 


LOAD  TYPE 


Load  Group 
1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


Primary 

1 

2 

3 

4 

5 

6 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 


Secondary  Coefficient 

K.  , 
i*l 

K.  _ 

1*2 

K.  . 

1,3 

K.  * 

1*4 

K.  c 
1*5 

K.  , 

1*6 


1 

2 

3 

4 

5 

6 
2 

3 

4 

5 

6 


Ki,7 

Ki,8 

Ki.9 

Ki,10 

Ki,U 


K 


K 


i * 12 
[i,13 
i * 14 


K. 

1*15 

K.  , , 
l , 16 


K 


i * 17 
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Load  Group 
18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


LOAD  TYPE  (Cont'd) 


Primary 

3 

3 

3 

3 

4 
4 

4 

5 

5 

6 


Secondary 

3 

4 

5 

6 

4 

5 

6 

5 

6 
6 


Coefficient 


Ki , 18 


Ki,19 
Ki , 20 


1,21 

Ki,22 

Ki , 23 

K.  14 
1,24 

K. 

l , 25 
K. 

l , 26 


1 , 27 


The  equations  for  the  primary  and  first-order  interaction  co- 
efficients, Equations  ( 29 ) — ( 34 ) , can  then  be  written  as. 


0£,i,n  Ai,£L£,£,n  + Bi,IL£,£,n 


where  for  £ * 1 - 6, 


i = 1 - 6 

I * 21  - (6  - £)  (7  ~ -)  + £ 
Ai,£  = Ki,£/Mi 


B.  T = K.  T/u. 
i,I  i,I  i 


(42) 


For  each  of  the  six  load  setups  Equation  (42)  can  be  fitted  for  the 

six  channels  and  n data  points  in  order  to  obtain  the  72  values  for 

A.  0 and  B.  For  the  six  sets  of  A.  „ and  B.  _ where  i = £ the 

i,£  i,I  i,£  i,I 

values  of  u . can  be  found  since  for  i = £,  K.  - = 1 and  u.  = 1/A.  . 

X 1 / ~ 1 1 f ^ 

The  first-order  interactions  can  then  be  calculated  from. 


K.  - = y . A.  - 

i, £ l l,  £ 


(43) 


\i  .B.  T 

l 1,1 


(44) 
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The  equations  for  the  second-order  interactions.  Equation  (41), 
can  be  expressed  as  follows. 


y . 0 


- K 


- K. 


- K. 


- K. 


K. 


i,£ 


l £,i,n  i,A  £,A,n  i,B  £,B,n  “i,C"£,A,n  "i,D"£,B,n 


**£ ,A,n^£ ,B,n 


(45) 


for  i = 1 - 6 and  where 


£ 


A B 


8-12  1 £ - 6 

14-17  2 £ - 11 

19-21  3 £ - 15 

23,24  4 £ - 18 

26  5 £ - 20 

C * 21  - (6  - A)  (7  ~ — ) + A 
D = 21  - (6  - B)  (7-^-B)  + B 

The  main  routines  used  in  the  balance  calibration  program  have 
been  listed  in  Appendix  B. 


MODEL  WEIGHT  EFFECTS 


Since  a combined  load,  second-order  calibration  is  used  to 
calibrate  the  balance  and  since  this  necessitates  referencing  all 
loads  to  a condition  of  zero  load-zero  balance  reading,  the  act  of 
just  placing  a model  on  the  balance  will  be  indicated  as  a load  in 
the  balance  readings.  With  a six-component  balance,  these  balance 
readings  can  be  used  to  determine  the  model  weight  and  center  of 
gravity.  These  values  could  then  be  used  to  calculate  the  balance 
load  caused  by  the  model  as  the  model  orientation  is  changed  in  the 
tunnel.  The  subtraction  of  these  loads  from  the  loads  measured 
during  the  run  would  then  leave  the  desired  aerodynamic  load. 

The  above  scheme  for  accounting  for  model  weight  effects  on  the 
balance  has  several  shortcomings.  First,  it  only  accounts  for 
static  effects.  Dependent  upon  the  model  mass  and  pitch  history, 
there  may  be  significant  dynamic  loading  effects  that  would  have  to 
be  corrected  for  either  experimentally  or  analytically.  Secondly,  a 
six-component  balance  must  be  used.  If  only  a four- component  balance 
is  used  and  the  model  center  of  gravity  is  not  on  the  balance  axes, 
it  is  not  possible  to  accurately  compensate  for  model  load  effects. 
Another  problem  that  arises  in  calculating  model  weight  effects  is 
that  the  balance  mass  also  affects  the  readings.  If  the  effects  of 
model  mass  on  the  balance  readings  are  to  be  calculated  as  a function 
of  model  orientation,  it  is  essential  (in  some  cases)  to  account  for 
both  the  model  mass  and  the  balance  mass  and  their  different  centers 
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of  gravity  independently.  Furthermore,  the  apparent  balance  mass 
that  influences  the  gage  readings  is  different  for  different  gage 
sections. 

Another  scheme  for  accounting  for  model  and  balance  mass  effects 
is  to  sweep  the  model  through  the  same  orientations  that  will  be 
used  in  the  data  run.  During  this  sweep,  record  the  balance  output 
channels,  reduce  these  data  to  loads  as  a function  of  model 
orientation,  and  then  subtract  these  measured  loads  from  the  data 
acquired  during  the  test.  The  disadvantage  in  this  technique  arises 
if  relatively  large  sting  and  balance  deflections  occur  as  a result 
of  the  aerodynamic  loads.  If  that  occurs  the  mass  effects  may  not 
be  known  for  the  same  exact  orientation  that  occurs  during  the  test. 

Neither  of  the  above  schemes  is  ideal.  Dependent  on  the  test 
parameters,  such  as  balance  mass/model  mass,  aerodynamic  load/model 
weight,  pitch  angles,  pitch  rates,  static  deflections,  etc.,  either 
of  the  techniques  may  be  preferred.  If  sting  deflections  are 
minimal,  the  technique  of  pitching  the  model  and  recording  the  mass 
effects  is  generally  preferred  at  the  White  Oak  Laboratory. 

TARE  MEASUREMENTS 

Tare  measurements  are  sample  recordings  of  all  test  parameters 
made  under  a no-flow  condition.  These  recordings  are  made  so  that 
a zero  reference  point  can  be  established  for  all  instrumentation 
just  prior  to  the  acquisition  of  data.  This  is  done  by  recording  a 
short  sample  of  the  balance  data  acquired  at  either  a fixed  pitch 
angle  or  during  a sweep.  These  readings  are  then  compared  to  the 
readings  acquired  during  the  weight  effects  sweep  for  the  same 
angular  orientation.  The  differences  between  the  balance  tare 
readings  and  the  weight  effect  readings  are  then  stored  as  balance 
offsets  which  are  applied  to  all  of  the  balance  data  recorded 
immediately  after  the  tare. 

If  any  change  in  the  balance  output  channels  occurs  due  to 
thermal  effects,  it  is  quite  important  to  obtain  a new  tare  in  order 
to  update  the  balance  reference  reading. 

DATA  REDUCTION 

In  the  data  reduction  each  line  of  data  (balance  data  recorded 
at  one  point  in  time)  is  reduced  independently  from  the  others. 

The  first  record  of  all  balance  channels  and  angular  positions  is 
read  from  the  data  tape.  These  data  are  then  adjusted  by  subtracting 
the  balance  offsets  that  were  determined  in  the  tare.  These  adjusted 
readings  are  then  reduced  to  loads  through  the  iterative  process 
described  in  the  section  on  General  Load  Resolution  Technique.  The 
aerodynamic  loads  are  then  found  by  first  interpolating  the  weight 
loads  in  order  to  get  the  proper  model  effect  at  the  same  orientation 
for  which  the  data  were  recorded  and  then  subtracting  the  model  load 
from  the  measured  load. 
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Once  the  loads  are  known,  it  is  possible  to  calculate  how  much 
the  balance  has  been  deflected  by  these  loads.  The  deflections  are 
calculated  from  the  information  determined  in  the  static  angular 
calibration  and  added  to  the  measured  orientation  in  order  to 
determine  the  true  orientation. 


These  loads  and  orientations  are  then  reduced  to  coefficient 
form  in  the  desired  axis  system  in  a separate  subroutine  called 
COEFF.  Generally  the  coefficients  are  defined  in  the  balance  axes 
as  follows: 


AFC  = C 

a 


YFC  = C 

y 


NFC 


RMC  = Cj_ 


Mx 

qaT 


PMC  = C 

m 


^Y_ 

QA£ 


YMC  =--  C 

n 


Mr 


QA£ 


where 

Q = dynamic  pressure  (psi) 

2 

A = reference  area  (inches  ) 

£ = reference  length  (inches) 

-Fx  = Fft  = axial  force  (lb) 

Fy  = yaw  force  (lb) 

-F„  = F„  = normal  force  (lb) 

L N 

= roll  moment  (in- lb) 

“y  = pitch  moment  (in-lb) 

M_  = yaw  moment  (in- lb) 
z 

The  various  axis  systems  and  transformations  are  shown  in  Appendix 
A and  the  data-reduction  subroutines  are  shown  in  Appendix  C. 
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GENERAL  COMMENTS 

The  purpose  of  this  report  is  to  explain  the  concepts  utilized 
at  the  White  Oak  Laboratory  in  making  static  force  and  moment 
measurements.  In  the  actual  utilization  of  these  techniques  many 
problems  arise  that  are  rarely  ever  brought  out  in  a general 
discussion,  and  unfortunately,  it  is  these  problems  and  how  they 
are  handled  that  determine  the  quality  of  the  data  obtained.  Some 
of  these  concerns,  such  as  the  accounting  ior  model  mass  effects 
and  balance  drifts,  have  been  alluded  to  already. 

In  general,  there  are  four  areas  that  must  be  examined  and 
checked  in  order  to  obtain  quality  data. 

The  first  question  is.  What  is  the  quality  of  flow?  Parameters 
of  concern  here  are  the  uniformity  of  the  flow  over  the  model  during 
the  test  period.  These  qualities  are  tunnel  dependent  and  outside 
the  scope  of  this  report  but  are  still  of  primary  importance  to  the 
value  of  the  data. 

The  second  area  of  concern  is  the  establishment  of  reference 
conditions.  Included  in  these  are  the  details  involved  in  determin- 
ing the  balance  output  with  no  aerodynamic  load  and  in  determining 
the  calibration  constants.  It  is  imperative  that  check  loads 
(multiple  combined  loads)  be  placed  on  the  model  and/or  balance  and 
reduced  in  order  to  check  the  validity  of  the  calibrations  and  the 
accuracy  of  pure  static  measurements  without  the  complications  of 
flow. 


The  third  area  of  concern  is  the  establishment  of  the  alignment 
of  the  model  and  balance  with  respect  to  the  flow.  Tnis  is  especially 
true  for  nonsymmetric  models  for  which  the  normal  force,  and  pitching 
moment  are  not  zero  at  zero  degree  angle  of  attack.  With  a 
symmetric  model  it  is  possible  to  move  the  model  in  the  flow  until 
the  orientation  is  determined  for  which  no  aerodynamic  force  or 
moment  exists.  This  location  establishes  the  point  where  the  total 
angle  of  attack  is  zero.  Depending  upon  the  accuracy  desired  and 
the  facility  in  which  the  test  is  being  conducted,  it  is  sometimes 
necessary  to  utilize  a symmetric  model  to  establish  these  conditions 
as  reference  conditions  for  a nonsymmetric  model.  For  models  which 
have  a slight  asymmetry,  it  is  possible  to  test  them  at  various  roll 
angles  (for  example,  0,  ±90,  180  degrees)  in  order  to  find  the 
orientation  of  the  model  wi *h  respect  to  the  flow. 

The  fourth  area  that  must  be  analyzed  carefully  is  not  always 
present  but  certainly  must  be  reckoned  with  in  some  small  hot 
facilities.  This  is  the  problem  of  how  to  handle  balance  drifts  and 
shifts  which  are  usually  caused  by  thermal  gradients  across  the 
balance  sections.  If  it  is  not  possible  to  insulate  or  cool  the 
balance  it  is  usually  necessary  to  limit  the  duration  in  which  the 
model  is  exposed  to  the  hot  flow.  This  may  mean  that  it  is  necessary 
to  obtain  data  for  just  a few  fixed  angles  rather  than  to  leave  the 
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model  in  the  flow  and  sweep  through  a large  range  of  angles.  If  a 
test  is  being  conducted  by  sweeping  the  model  through  large  angles  of 
attack,  the  thermal  effects  are  angle  dependent  and  may  not  be 
evident  to  the  project  engineer  unless  multiple  sweeps  are  made  and 
the  data  are  compared.  It  is  possible  for  the  thermal  effects  to 
follow  the  angle  of  attack  closely  enough  so  that  the  differences 
between  tares  conducted  before  and  after  the  run  are  very  small  but 
yet  the  effect  on  the  data  at  high  angles  may  be  quite  large. 

The  ideal  way  to  eliminate  the  above-mentioned  problems  would  be 
to  improve  the  design  of  the  balances  and  tunnels,  but  since  this  is 
not  always  financially  possible,  the  project  engineer  must  resort 
to  procedural  changes  in  the  manner  in  which  the  test  is  conducted. 
This  may  mean  that  additional  flow  alignment  runs  are  necessary,  or 
that  sweeps  may  not  be  possible,  or  that  additional  tares  may  have 
to  be  taken.  All  of  these  things  must  be  considered  in  the  initial 
design  of  the  test  program. 
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APPENDIX  A 
AXIS  TRANSFORMATIONS 

Within  the  daca-reduction  process  there  are  a number  of  axis 
systems  which  are  used.  The  tunnel  axes  are  considered  to  be  the 
reference  axes.  These  axes  may  be  defined  slightly  different  in 
each  of  the  tunnels  but  in  general  both  the  XT  and  axes  are 

defined  to  be  in  a horizontal  plane  (established  with  a precision 
level)  with  the  XT  axis  pointing  into  the  flow  and  the  vertical  axis, 

ZT,  pointing  down.  All  angular  rotations  of  the  sting  are  defined 

with  respect  to  the  tunnel  axes  through  the  angles  0,  tp,  and  <p, 
and  generally,  in  that  order.  Since  the  tunnel  axes  are  somewhat 
arbitrary,  the  matter  of  real  importance  is  the  location  of  the 
sting  with  respect  to  the  wind  vector.  The  small  angles  A0W  and 

At/>w  are  used  to  define  the  tunnel  axes  with  respect  to  the  wind  axes 

(See  Fig.  A-l) . 

The  transformation  matrix  to  go  from  the  wind  axes  to  the 
tunnel  axes  in  the  order  of  pitch,  yaw,  and  roll  is  called 

cosA^wcosA0w  sinAi|>w  -cos  A^wsinA0w 

sinAi^wcosA0w  cosAi^w  sin  A4>wsin  A0W 

sinA0  0 cosA0 

w 

where  the  order  is  A0W  -*•  A^w  A4>w . 

The  transformation  matrix  to  go  from  the  tunnel  axes  to  the 
sting  axes  is  called  [£.T£,]. 

cos^cos© 

[£.  1 = -cos<j>sinii/cos0  + sin4>sin0 

X b 

sin<f>sinij>cos0  + cos$sin0 

where  the  order  is  0 -*■  i|>  -*•  $. 

If  a dogleg  sting  is  used  or  if  the  first  rotation  is  in  the  yaw 
plane,  a different  transformation  matrix  must  be  used  to  go  from  the 
tunnel  axes  to  the  sting  axes.  The  sting  orientation  would  be  as 
shown  in  Figure  A- 2.  The  transformation  matrix  would  then  be. 


sintj'  -cost^sin© 
cos4>cos^  cos4>sin^sin0  + sin<J>cos6 
-sin4>cosi^  -sin4>sin^sin0  + cos<f>cos0 
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[£tsy] 


■ COS0COSl|< 

-cos4>sini|<  + sin<J>sin0cost|; 
sin<j>simj;  + cos4>sin6cost|/ 

where  the  order  is  ip  -*■  0 -►  <J>. 


cosQsini^  -sin0 

cosQcosip  + siwpsindsiwp  sin$cos0 

-sin<j>cos<J;  + cos<{>sin0sini|'  cos0cos<)> 


The  variation  in  the  balance  orientation  with  respect  to  the 
sting  is  due  to  the  deflection  of  the  balance  and  sting  due  to  the 
loads.  These  deflections  are  defined  as  6,  y,  and  e (see  Fig.  A-3) 
The  transformation  matrix  is 


"sn1 


I'SB1' 


cosycosS  siny 

-cosesinycos6  + sinesin6  cosecosy 

L sinesinyccs6  + cosesin6  -sinecosy 


-cosysinfi 

cosesinysin6  + sinecosfi 
-sinesinysinfi  + cosecos6j 


where  the  order  is  6 + y ■»  e. 

If  the  model  is  rolled  with  respect  to  the  balance  through  an 
angle,  <J>„,  the  transformation  matrix  is  [£___]  (see  Fig.  A-3). 


<W  - 


o 

COS<j> 

-sin<J> 


0 

sin<f> 

cos<j> 


When  the  model  moment  reference  center  is  not  located  at  the 
balance  center,  it  is  necessary  to  translate  the  balance  moments. 
The  distances  along  the  rotated  axes  are  called  DXMC,  DYMC,  and 
DZMC  and  are  positive  if  the  moment  reference  center  is  located 
forward,  right,  or  down  from  the  balance  center  (see  Fig.  A-4) . 
The  moments  referenced  to  the  model  moment  reference  center  are 
calculated  as  follows. 


V 

’dxmc 

“yc 

X 

«Y' 

+ 

DYMC 

.V 

DZMC 

N 


-pi 
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FIG.  A4  RELATIONSHIP  BETWEEN  BALANCE  CENTER  AND  MODEL  MOMENT 
REFERENCE  CENTER 
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APPENDIX  B 


BALANCE  CALIBRATION  PROGRAM 
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BALANCE  CALIBRATION  PROGRAM 
APPENDIX  B 

PROGRAM  BALCAL  (lNPUT=F4*0DTFUr  = b4  »TAPE.5=INPUT  tT4PEb=0UTPUT  ) 


THIS  PROGRAM  IS  USED  TO  REDUCE  HALAPCc  C AL  IHhAT  I OK  LOADINGS 

THE  FOLLOWING  INPUTS  ARE  DEFINED  AS 

HAL  = RAPE  OF  BALANCE 

PC  = MMHFR  CF  HALANCS  CCN-hONEMS 

ITYPE  = 1 INDICATES  A POPEM  PALANO 

= ? INDICATES  A ECHCE-PCmNT  FiALANCt 

OATE  = OATF  CAIIHHAIIOF  t*AS  PtRFCwPED 
F.NG  * ENGINEER  PERFORMING  CALIEmaUON 
OAHE  = DATA  ACQUISITION  SYSTEP  UStL' 

SETSCAL  = NOP  INAL  mV  FOI  L SCaLF  SL N bl T 1 V I T Y . . . .NOTt  IF  THIS 
IS  A REDUCTION  OF  A MO«F.NT  HALAULt * ITYPE  = 1.  THE  SETSCAL 
VALUES  FOR  CPANNFL  1 MUST  EQUAL  THE  VALUES  FOR  CHANNEL  2 
AND  CHANNEL  2 RUST  PATCH  ChANNtL  A. 

OFF  * INITIAL  REAPINGS  OF  ALL  CHANNELS  *IlH  NO  IOAI) 

SCALE  = SENSITIVITY  SETTINGS  Af  *HlCH  l OAUb  HERE  PEASUHEU 
*T  = VALUES  CF  LOAUS  HUNG  (Lh  Ch  IN-Lb) 
h = DARE  CHANNEL  READINGS  (CULMS) 

L = LO AO  TYPE  NUMHEh 
P » CHANNEL  NUPHER 
N a LOAD  NL'MEEFi 

LS  = 11  AT  FND  OF  1NCIVIIUAL  LCADIMj  bETLP 
* 95  AT  _NU  OF  CAL1HHAT10N 


DIMENSION  WT(27.b»l<)>.4<b.27>»H<b»27).DC(b.27)»X(10),Y<iU>»AN<2)* 
1PRMSEN  <b)  .OFF  (b) .HD (2  7 «b  « 1 0 ) .SCALE  (fc*?7)»Fi(27«G»10)tT11LE  ( 27 ) » 

1FL  ( 6 > .PL  (36)  • C 1L  (fc>  .C?P  (ft)  .ThLN  (b)  .LELT  (b)  .NN(?7)  .SFTSCAL  (6) 
DIMENSION  RAL  (2)  .CATE  (?)  .FN(-  (?) 

HEAL  NH(b) 

DATA  TITLE/?HFN.2HFY.2HMY,2Hi,/.2HMX«2nH  A.4HFnFN»4HFNF Y.4HFNPY. 
14HFNM7.4HFNMX.4FFNFA.4HFYPY.4hFYMY.4HF YMZ.4RE YMX »4HF YF A ♦ 4HP YM Y . 
14HMYPZ«4HP  YMX  .4FPYF  A.4hp7M7 ,4Hm7MX  .4ht"/F  A » 4HM XM X . 4hp XF" A »4HF  AF  A/ 


* ARRAYS  ARE  INITIALIZED 


Of)  P P = 1 » F 

P OFF (P)=THCK (P)=DEL T CM) aPRMSE N ( m) *0 . U 
DC  9 L=1 *27 
NML  1=0 
DO  V P=1 »E 
A (M.L ) =0.0 
SCALE (M.L ) =0.0 
0C  (M.L) =0.0 
DO  9 N = 1 » 1 0 
9 PT  (L»P*N)=0. 

• 

* all  CALIHRATION  DATA  AWE  HEAD  in  Ht  HE 


10  READ  05.8  0 0?)  baL.mC.ITYhE 
REAO  05.8  0 0 2)  DATE 
RE  AO (S  « BOO? ) ENG 
RF  AD (5*800?)  DARt 

READ  (S»800S)  CUP.rt'M.  (SETSCAL  O')  .»  = ) .b) 
800?  FORMAT  (2Alrt.TS.T5) 

20  RF  AIKS. MOOS)  I.  .DU*  . (OF F <P  ) .P  = J .h ) 

30  READ (S .HOOS ) l.  DIM. (SCALE (P .1 > »P= 1 »b> 
DO  35  Mal.MC 

35  OFF  (P  ) =CEF  <P)  «SC  ALF  (M.L  ) /SF.  1 SC  «L  (P  ) 

40  RFAO  (8.80CS)  L.N.  (*T(I  .P.M  .M=l.b>  .LS 
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RAL00100 


8AL00110 


RALUOlbU 
BALOOlbO 
RAL001 70 
R AL  00 1 80 
HAL00190 
HAL00200 


HAL002JO 
RAL0024U 
H ALo02bO 
HALO0?bO 

PAL  002  7 0 
P AL002M0 
MAL00290 
PALOOJOO 
8ALOO310 
PAL00J2U 

PAL00340 

PALOOJbO 


HAL00390 
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APPENDIX  B (CONT.) 

0005  FORMAT  (215. 6F10. 0.12)  PALU0410 

• PAL0042U 

• LS  = 11  AT  END  OF  InCIVICUAL  LC  AL)  I NO  btlLPS 

• LS  * OS  AT  END  OF  CALIF*  CATA  PAL00440 

• PALU0470 

50  READ (5  *0005)  L »N  * (h  < L »M  »N ) . M = 1 . 6 ) »LS  HALOU40O 


* READINGS  ARE  NO*  ACJLSTED  FOR  PROPER  SCALE  SETTINGS 

« 

52  00  55  Nsl.MC 
NH (M) *H (L  »M»N) 

NH  (M) =NH (N ) *SCALE  <P .L) /SE 1 SC aL  CM) 

55  CONTINUE 

* PAL00490 

* MOMENT  GAGE  READINGS  A RE  ADJUSTED  FOR  UTILISATION  AS  FPROE  HALuOSOO 

* AND  MOMENT  GAGES  IF  NECESSARY,  AM)  ADJUSTMENTS  ARF  MADE  FOR  OFFSETS 

* HAL00S20 


IF (ITYPE .FO.l)  CO  TO  60 
GO  TO  02 


60 

RO  (L.l.N)  =NH(2)-0FF  (2)-  IKE-  ( I > -OFF  ( 1 ) ) 

RD ( L « 2 »N  ) *NH ( A ) -OF  E (4)-(NP (2) -OFF (3) > 

HO  CL»3»N)*N'H<  1 )-OFF  (1 ) *NH  ( 2 ) -OFF  ( 2 ) 

RD(L»4.N)*NH(3)-0FF(3)  ♦NRU)-0EF(4> 

RD(L.5.N)*NH(5)-0FF<5) 

RD (L  *6  »N ) *NH (6  > -OFF  (6) 

GO  TC  70 

62 

DO  65  M*1.MC 

65 

RD (L»m,M*nh<m)-0FF ( M ) 

70 

NN ( L ) =N 

PAL  U05V0 

IF (LS.tt.ll)  GC  10  2o 

RAL00S94 

IF (LS.F0.S9)  GC  TO  PO 

PALOObVb 

GO  TC  40 

MALUObOO 

• 

PALOObl 0 

* 

SOLUTION  FOR  PRIMARY  SENSITIVITIES.  E IRsT  ORDER  JMERACTIUNS. 

PAL  00620 

• 

ANr  SECOND  OROFP  LCAC  SOD  ARE  I’  INTERACTIONS 

PAL00630 

ft 

PAL006AO 

PO 

CONT INUE 

PALoOhSO 

00  TC  400 

PAL  00652 

*0 

CONTINUE 

PAL00634 

DO  200  L=1»MC 

0 AL  00660 

00  ISO  1*1. RC 

PALO0670 

NO=NN (L) 

PALU06HO 

IE (NO.EQ.O)  GO  TC  200 

PALO0690 

DO  150  N* 1 .NO 

PAL00700 

X(N)=4T(L,L»N) 

9ALUU710 

YIN) =R0 (L . I .N ) 

PALOO  720 

150 

CCNTINUt 

HAL00730 

ft 

PAL00740 

ft 

ECUATICN  Y=A < 1 >*X*A (?) «x**2 

PAL0O750 

ft 

PALOO  760 

CALL  F1T1  (NO.x.Y.AN) 

PALO0770 

ft 

PAL007B0 

A ( I • L ) =AN  ( l ) 

H ALuO  790 

I 1=21- ( ( 6-L ) * ( 7-L ) /? ) *1 

RALOOfcOO 

0(1. II)* AN (2) 

PALOOfc 1 0 

190 

CONTINUE 

PAL  00020 

200 

CONTINUE 

PAL00630 

ft 

PALOO04O 

ft 

PRIMARY  SENSITIVITIES 

RALO0M50 

ft 

° AL  00060 

00  210  1*1, PC 

PAL008/0 

8C ( I » I ) = 1 • 

PAL00680 
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210 


It  (A < I » I ) .EQ.O.C)  00 
PRMSEMI)  = 1./A(I,I) 
CONTINUE 


TC 


APPENDIX  B (CONT.) 

210 


227 


229 

230 


FIRST  CRDFR  INTERACIIONS 

00  230  L=1»MC 
DO  229  1=1, MC 
IE(I.EO.L)  GO  TC  227 
BC  ( I .L  ) =RRMSFN  < I ) * A < I .L  ) 

SECONO  ORDER  LCAO  S'lLAPFC  INTER  ACTIONS 

11=21- ( (6-L)«(7-L)/2)*L 
BC ( 1 ♦ 1 1 ) =PRMSFN ( I ) *B ( I » 1 1 ) 

CONTINUE 

CONTINUE 


* SOLUTION  OF  REMAINING  SECOND  UNDER  INTERACTIONS 

DO  300  L=P,2E 
NC=NN (L  ) 

IF(NO.EO.O)  GO  TC  300 
IF (L.GE.B.AND.L.LF.l?)  23c.237 
235  L A = 1 
LR=L-6 
GO  TC  270 

237  IF (L.E0.13)  300.241 

241  IF<L.C-fc.l4.AN0.L.LF.17>  243,245 

243  LA=2 

LB=L-  1 1 
GO  TC  270 

245  IF(L.KO.IC)  300,247 

247  IFCL.Gf .19.ANn.L.LF.21>  249,251 

249  L A=3 

LB=L-15 
GO  TO  270 

251  IF (L.FU.22)  300.253 

253  IF  (L.C-E  .22.AN0.L.LF.24)  255.257 

255  LA=4 

LR=L-ld 
GO  TC  270 

257  IF(L.E(1.25)  300,259 
259  LA*5 

LH=L-20 


270 


LC=21-(  (6-LA)*  17-LA  )/2)  *1.4 
LD=2 I - ( (b-LB)«(7-LB)/2)*LH 
DO  290  1=1. MC 
DO  2B0  N= 1 .NO 

Y (N)  =P9MSFN(  I ) •RfHL.I  .O-HC  ( I , l.A ) **T  (L.LA.M  -BC  ( 1 »Lb)  #*T  ( L , Lb  , N ) 
1 BC  ( I »LC)««T  (L.L  A.N)**?-ML  ( I » L l> ) * h I (L.LH.N)  **2 
280  X(N)=»T(L»LA»N>«i»T  (L»l  h.N) 

CALL  F IT2 (NO.X.Y .AN) 


Y=A (1 )*x 


ECUATICN 

290  BC ( I »L ) = AN ( 1 ) 
300  CONTINUE 


* CALIBRATION  LOADS  ANC  READINGS  11STEC 
GO  TC  552 


BAL00B9U 
PAL 00900 
PALU0910 
BAL00920 
B AL00930 
BAL0U94U 
RAL00950 
BAL00960 
BAL00970 
BAL009B0 
BAL00990 
BAL01000 
FAL01010 
BAL01U20 
HAL01030 
ft ALO 1 040 
HALO 1 050 
BAL01060 
PALO  1 0 70 
BALOlOoO 
BAL01090 
BALul 1 00 
HALO  1 1 1 0 
PALol IcO 
BAL01130 
PAL01 140 
BAL01150 
PALUl ItoO 
ftALUll  70 
RALOl 1B0 
B AL  0 1190 
HALO  1200 
HAL01210 
HAL01220 
HALO  1 230 
BALul240 
ft AL  0 1 250 
PAL012G0 
PAL01270 
BALol2e0 
PALO  1 290 
PALol JOO 
PALU1310 
PALUl 320 
PAL0I330 
PAL0134U 
RALOl 350 
ft ALO 1 3ft0 
BALOI37U 
PALo 13dO 
*B ALO 1 390 
HALO14U0 
PAL014I0 
BAL01420 
PALo 1430 
PALol 44  0 

palo i 450 

PAL01450 
PAL01470 
HALO  1 4B0 
HALO  1 490 
BAL01500 
BALUl 502 
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APFENOIX  B (OONT.) 


400  WHITE  (6.9001 ) 

9001  FOF<HA  I (1H1»53X*NAVAL  ^UPFAO  wFAHCNS  CENT E P*//5 ?X*wh IT t OAK 
1T0HY* ) 

WHITE  (6,9002)  HAL. ENG 

900?  FOFJHAT  (/4X*HALANCF  *?  A 1 0 . Ml  A*t  *C- 1 NE  E H *2A1U) 

■HITE (6.9003)  OA IF »UAWF 

9003  FOHHAT (4X«CAL IHfiAT ION  CATE  *?A 1 0 • b6X*UAKE  **A10> 

■HITE (6.9040)  SETSCAL 

9040  FCHHAT (/52X*HV  FULL  SCALE  SE  TT  J NGS  *fc  IF  10.2)/) 

WHITE (6.9004) (TITLE (L).L=1.6)»(TITLk ( L ) »L  = 1 .6) 

9004  FOHHAT  </*  L N*6 ( 8X A? ) »6X»6(PXA2)  ) 

■HITE (6.9005) 

9006  FOHHAT  ( 1 0X.6  (6X«LCA0*)  .6 A. 6 (7X*HD(:*)  ) 

00  500  L=1.27 
NO=NN (L ) 

IF (NO.EO.O)  GO  TC  600 
00  499  N=1.N0 

■ H IT  E ( 6 « 90  06 ) L.N.  (*T  (L  ,M.N)  .H  = J .6  ) . (MML.H.N)  »H*1.H) 

9006  FCHHAT (/2I4.4X.6F1 0.4.3X.HF 10.1 ) 

499  CONTINUE 

500  CONTINUE 
■HITE (6.9001 ) 

■ PI  IE (6.900?)  HAL. ENG 
■HITE (6.9003)  CATE.UAhE 

■ wITE  (6.9007)  (TITLE  (L)  »LS  1 .6  ) ♦ (L.L  = 1 .6 ) 

9007  FOHHATI/*  L N«6 (H*A?)  .6X.6 (HX »H*  1 1 ) ) 

DO  550  L=l,?7 

NC=NN (L) 

IF (NC.FU.O)  GO  It  660 

■HITE (6.9040)  (SCAIE(h.L).  H=1.HC) 

00  549  N=1,N0 

■HITE (6,9006)  l.N.(wT(L.H.N),H=l»6)»(Fl(L.6,N)»H=l,6) 

549  CONTINUE 

550  CONTINUE 

GO  TC  90 
65?  CONTINUE 


PALOlblO 
PAL0152U 
LAbOHAHALOlbJO 
HALO  1 54  0 
HALU1560 
HALolSbO 


PAL01690 
HAL0160U 
PAL0161U 
HAL  0 1 620 
HAL01630 
HALO  1 64  0 
HALO  1650 
HALO  1660 
HAL01670 
HALO  1660 
HAL 01690 
HALOl 700 
PALO  1710 


H AL  0 1 7**  0 
HAL 01750 
HAL*1 1760 
HALOl  1 10 
PAL01760 

HALOl  / 9 0 
HALOlbOO 
bALUlblO 
PALO  1 H?0 
HALOl  t)3o 
PALolbJ7 
PAL  0 1 M Jh 


* CALIPKATION  CONST  AN  I b llSTEU 


HAL01H40 
HALO  165  0 


■HITE (6.9001) 

■ HITE  (6.900?)  HAL  ,FN(- 
■hITE (6.9003)  CaTF.Dahl 

■ HITE  (6.901'H)  HALOlbOO 

9CCP  FCHHAT  ( / ,13X*CAl  THHATIO  SYSTFN  V ALUF  S*  JOA*nOHNA|  TZtU  1MEHACIIC 

1 N */ 74 x ‘COEFFICIENTS*) 

■H I T t (6,6009) (NK ,NK=1 .6 ; HALOl 690 


9009  FCHHAT  (/I  1X*HCHENT*4X*CAH6  it  At.F  *e6*.6I  12) 

■ HITE  (6. 9010)  (TITLE  (Nk  ) »NF  = 1 ,t> ) 

9010  FCHHAT  ( i?X*GAGF«6X* (HV-FUI  L SC ALE ) «2?A ,6 ( 1 OX A? ) / ) 

UC  310  NK=1,6 

■ hITE (6.9011)  NK.SFTSCAL (Nr ) . r N . TI  ILL  (NK) » (HC(I.NN) .1  = 1.6) 

9011  FCHHAT  (14X,I1.3XcH.?,i/X.Ill,iXA?.jA,cti2.5) 

310  CONTINUE 

DO  315  NK=7,in 

■PITE(6.9017)NK.TITLE(NK),(6C(I.NK),I=l,6) 

9017  FOHHAT  < 52X  , I?  • 1 X . A*V  , 1 X . 6F  ]f  .6) 

315  CONTINUE 

■HITE (b»90?l ) (FC(I.11)»I=1.6) 

90?  1 F0HHAT(1?X.  «PwlHt  SENS  I T I V I I 1 1 5*21 A*i 1 F NHX*  l X . At  1 2 . 6 ) 

■ H IT  F (6.90??)  (Ff.  ( I » i?)  » 1=  1 .6  ! 


PAl.  01910 

PAL01 9J0 
HA 

MAL0196U 
PAL 0 1 9 i‘  II 
BAL019H0 
HAL  0 1 990 
PAL02000 
BALO2010 

PAL0?030 
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9022 

APPENDIX  B (CONT.) 

FORNATUOXMLP  CH  IN-LH  PER  COUNT ) *19X*12  FNFA*1X.6F12.5> 

WHITE (6.90P3)  <PC(I. 13) .1  = 1.6) 

PALU2050 

9023 

FORHAT (52X*13  FVFV*1X,6F 12.5 ) 

PAL0206U 

DO  330  NK  = 1*6 

PAL0P070 

NZ«NK*13 

H AL02080 

WHITE (6.9024)  TITLE <NK) .PHkSEN (NK).NZ.IITLF (NZ) . (HC(I.NZ) 

,1=1,6)  RAL02090 

9024 

FORHAT (12X.A2.5X.E 12.5 . 2 IX • 12. 1 XA4 * 1X.6E12.5 ) 

PAL021 00 

330 

CONTINUE 

halo?: 10 

DO  340  NK  = 20 .27 

PAL02120 

WHITE (6.9030)  NK.T ITLF (NK ) . (HC ( 1 .NK ) . 1 =1 .6 ) 

PAL  021 30 

9030 

FOPHAT  (52X.I2.1X.A4. IX, t;E12.S) 

8AL02140 

340 

CONTINUE 

BAL02150 

* calibhation  CHECK 

BAL0?16n 
fi AL  02 1 70 

WHITE (6.9001) 

WHITE (6,9002)  PAL.FNP 
WHITE (6,9003)  CATE, DANE 
WHITE (6. 9038) 

BAL02180 
BAL021 90 

9038 

FORHAT ( / ,60  X*C AL IHHAT ION  CHECK*// 

B AL02200 

134X*APPLIFH  WEieHTS*34X«JlFFEHENCE  BETWEEN  LOAD  HEADING  AND  CALCUL 

2ATE0«/84X«RE ACINU.  K(6X27)  X L(27xl)  IN  COUNTS*) 

9041 

WHITE (6.9041 ) SETSCAL 

FORK AT (/59X* (HV-FLLL  SCALE ) *6 (F 1 0 .2 ) ) 

WHITE  (t.,90  04)  (TITLF  (L ) »L=1.6)  . (TITLE  (L>  »L*1.6) 

PAL02230 

WRITE (6.9005) 

PAL022h0 

DO  1000  L= 1 .27 

HAL02250 

N(’  = NN  (L  ) 

HAL02260 

IE  INC. EG. 0)  00  TC  1000 

HAL02270 

DO  999  N=l.NO 

PAL02280 

LL  =6 

BAL0229U 

DO  530  1=1.6 

PAL02300 

530 

FI  (I)=wT (L.I.N) 

PAL0P310 

DO  540  1=1,6 

RAL02315 

DO  540  J=  I , 6 

PAL02320 

LL-LL* 1 

PAL023J0 

540 

PL(LL)=FL(I)*FL(J) 

HAL02340 

DC  560  1=1.6 

PAL  02350 

C?P(I)=0. 

PAL02355 

560 

C 1L  ( I ) =0 . 

PAL02360 

DO  600  J= 1.6 

8 AL02370 

DC  600  1=1,6 

PAL02380 

600 

C1L(J)=BC(J,I)*FL(I)*C1L(J) 

BAL02390 

DO  700  J=  1.6 

PAL02400 

DC  700  1=7.27 

PAL  0?4 1 0 

700 

C?P  (U)  =C2F  ( J)  + HC  U » I >*HL  ( I ) 

PAL  02420 

DO  800  J=1 .6 

PALOP4J0 

IF (PHHSEN ( J)  .FQ.O.O)  GC  TC  «00 

HAL02440 

THCK ( J ) = ( C ) L (w) «C?P ( J) ) /PH^SEN (U) 

HAL02450 

DELT (J) *HC (L .U.N )-THCK (w> 

PAL02460 

800 

CONTINUE 

PAL  024  7 0 

WHITE  (6,9006)  L.N.  (WT  (L  »H»N)  »C  = 1.6)  . (L'ELT  ( J)  » J=l»b) 

BAL02480 

999 

CONTINUF 

PAL0P4V0 

1000 

CONTINUE 

HAL02500 

END 

PAL02510 

SUBHCL  T INF  FIT1  (NPT.XX.Y.Ar, ) 

PAL02620 

01  PENSION  XX  (10)  »Y  ( 10)  . AN  ( 2 ) *I»(S)  »A  ('-»*  iO>  ,A(26,3) 

HAL02530 

EXTERNAL  VFl 

PALO2540 

N < 1 ) *2 

PAL0P550 

N (2) =NHT 

PAL0P560 

N (3) *2 

PAL025  / 0 

N (4 ) =2 

HALO2560 

B-fi 


i 
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N<5)=0 

PALU2590 

00  10  1 = 1.2 

BAL026U0 

A ( I * 3 ) = I 

HAL02610 

10  A ( 1 1 1 ) =0 • 

HAL02620 

PAL02630 

00  20  L=1.NPT 

BAL0264U 

X ( 1 *L ) =XX (L  > 

PALU2650 

X ( 2 »L ) = Y (L) 

HAL02660 

20  X (3.L ) = 1 • 

HAL02670 
P AU  02680 

DFLT  = 1 • 0E-6 

HAL02690 

CALL  LSQSLRIN.X.rt.VFl.OELT) 

HAL02700 

00  30  1=1.2 

P ALU?  7 1 0 

AN (I)sA(I.l) 

PALU2720 

30  CONTINUE 

PALU2730 

RETURN 

BALU2740 

END 

PALU27S0 

SUBROUTINE  VF1 (N  *X .DF  * A ) 

PAL  02760 

DIMENSION  OF ( 1 ) . X ( 1 ) • A ( 1 ) ,N(1) 

PAL 027  7 0 
PAL  02760 

EOLATION  Y = A ( 1 ) A. A (?) X**2 

PALU2790 
P AL02600 

OF  <1)=X(1) 

PAL02H10 

OF  (2)  =X  ( 1 ) **? 

PALU2820 

X (9) =A > 1 ) • )F ( 1 ) *A (2) *0F (?) 

BAL02H30 

RETURN 

RAI.02B4U 

FND 

PALU2MS0 

SUBROUTINE  FIT2INFT.XX.Y.AN) 

P A l U2M60 

DIMENSION  XX ( 1 0 ) *Y (10) «AN(1)  *N  ( 5 ) * X ( V • 1 0 ) .A126.3) 

PAL02M70 

EXTERNAL  VE2 

PAL02680 

N ( 1 ) = 1 

eAL0?P90 

N ( 2 ) =NPT 

HALU29U0 

N ( 3 ) =2 

PAL02910 

N ( 4 ) =2 

HALU2920 

N (S) =0 

PAL0293U 

A ( 1 .3) =1  . 

PAL  02940 

10  A(1.1)=0.0 

P AL029SO 

DO  20  L=1.NPT 

PAL  029c  0 

X(l,L)=XX(l) 

PAL  0297  0 

X ( 2 . L ) = Y ( L ) 

P ALU2960 

20  X ( 3 . L ) = 1 - 

PAL029VU 

UtLT=1.0E-6 

R AL  U 3000 

CALL  LSGSLh(N.X,A.VF2.0ELT) 

PALU3010 

AN ( 1 ) =A ( 1 • 1 > 

PAL  03020 

RE  TURN 

BAL03030 

END 

PAL  03040 

SUBROUTINE  VE 2 (N . X . OF ♦ A ) 

RALU30S0 

DIMENSION  OF ( 1 > .X( 1 ) .A ( 1 ) .N ( 1 ) 

PAL03060 
PAL  030  70 

EOLATION  Y = AX 

PAL03080 
HALO 30V0 

of  < i >*x  m 

PAL03100 

X (9) =A ( 1 ) •OF ( 1 ) 

RAL031 10 

HF  TURN 

RAL03120 

END 

PAL03130 
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DATA  REDUCTION  SUBROUTINES 
APPENDIX  C 

SUBROUTINE  BALPP  (JCMPT5.HP.C11.C1  LC2 ) FOhCE 

PURPOSE 

CALCULATE  THE  INVERSE  NAIHIX  OF  THfc  NORMALISED  1ST 
ORDER  INTERACTIONS  *N0  TFF  HhCCUCI  OF  THIS  INVFRSE 
MATRIX  * I TH  THE  NORMALIZED  2ND  ORDER  INTERACTION  MAIHIX 

DESCRIPTION 

JCPRTS  INPUT  NUMBER  CF  RAlANCE  COMPONENTS 
HP  INPUT  NChmai  I2EO  PALAN’CE.  PAIhIX 

ClI  OUTPUT  INVERSE  R A TH 1 A CF  1ST  ORDER  I NT  F HAT  T 1 UNS 

C1IC2  OUTPUT  PRODUCT  OF  ClI  ANo  2ND  OHOEH  INTERACTIONS 


DIMENSION  C2  ( lift ) .H  <b.  1)  •(  1 IbiM  .HP (N.2Y)  .Cl  I OH)  .Cl  IC£(lt-2> 
M*JCPPTS 
NCsM«<M*1)/2 
00  30  J*1.M 

B( Jil)*1.0 
DO  30  Ks]«M 
30  Cl (J.K)*HP ( J.K) 

Call  matrix (ci .p »f .d.o.ceterm. in ) 

L = 0 

00  AO  K* 1 • M 
DO  40  J* 1 * M 
L*L*1 

40  ClI (L ) =C1 (J.K) 

L = 0 

K 1 *P* 1 
K?*P*NC 

DO  50  K=K1.K2 

DO  50  J=1.M 
L=L*  1 

50  C?(L)=BPU.K) 

#*# 

CALL  (jMPXFAtClI.CE.Cl  IC2.M.P.NC) 

RF  TURN 
t NO 


FOHCt 

FOHCE 

FOHCt 

FOhCE 

FOHCE 

FOHCt 

FOhCE 

FOhCE 

FOHCE 

FORCE 

FOHCE 

FOHCt 

FOHCt 

FOHCt 

FOHCE 

FOHCt 

FOHCt 

FOHCE 

FOHCt 

FOHCt 


FOHCE 

FOHCt 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  I TEH A ( TH »LO AO .HHPSF N .C II .C 1 IC2 ♦ JCmrTS  . I I P A A » ACCUh . 1 FH ) FORCE 
PURPOSE 

REDUCES  THE  CATA  TO  LOADS  PY  USING  ThL  BALANCE  MATRIX  AND 
ITERATING 

DESCRIPTION  OF  PARAMETERS 

TH  CFiANNEL  READINGS  IN  COUNTS 

LOAD  LOAD  COPPONFNTS  CORRECTED  For  INTERACTIONS 
(FOPCF  AND  MOMENT  UMTS) 

ClI  INVERSt  OF  N0HPALI7F0  1ST  OROF  H INTtHACTIONS 
C 1 ICE  PRODUCT  OF  rn  anC  NOHPALI/tD  2ND  UHOFR  IMtHACTIONS 
JCMRTS  NUMBER  CF  FiALANCF  COMPONENTS 
1TPAX  MAX  I ML M NUMHBH  CF  TTFHATIOnS 
ACCUR  SMALLEST  LOAD  UNIT  T BAT  CAN  Ut  wESOLVtU 
IEH  ERROR  INDICATOR  FOR  CON  VtM->t  NCt 
I EH  * 0 LOADS  CONVERGED 

IEH  = 1 LOADS  (>IC  NOT  CUNvF  HGt  HEFUHE  I TPAX 
HAS  REACHED 


C ~'l 


t 


o r> 
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APPENDIX  C (CONT.) 

HE  AL  LU<6> .PWMSEN (ft) .TH (ft) .EHSO(ft  > . H ( 3ft  ) . C 1 I <36  ) . Cl  I C2  ( 1 ftc ) .LC  Al)  (6  FOHCE 


1) .DELTA (ft) . ACCUF (ft) .E PS  I (ft)  FOHCE 

M=JCMPTS  FOhCE 

N*M* (M* 1 ) /?  FOHCE 

00  20  1*1 « JCMPTS  FOhCE 

20  LO(I)=PRMSFN< I)«TF  (I)  FOhCE 

• ••• 

CALL  6MPXF  A (C  ] I.LL.LOAD.M.m,  1) 

L=0  FOHCl 

00  ftO  1*1. M FOhCE 

EPSO ( I ) *0 • 0 FOhCE 

DO  60  J*  I » M FOHCE 

L=L* 1 FOHCt 


IF(I.NE.J)  GO  TC  59 
IF  ( L 0 A 0 ( I ) ) S».5S».59 
5d  P(L)=-LOAC(I)«LCAC(I) 
GO  TO  bO 

59  P (L)  *LOAD  ( I ) *LOA()  ( J) 

60  CONTINUE 


I T*0  FOHCE 

DO  300  M*  1 . I TMA X FOHCE 

IT*IT*1  FOHCE 

80  CALL  GMPXFA (Cl IC2.H.FPSI.M.N.1.  ) 

90  CALL  GMS8F < EPS  I » EPSO.OF L T A .M ) FOHCE 

100  CALL  GMSHF (LOAD. DELTA. LOAD, ft)  FOHCE 

L=0  FOhCE 

00  105  1*1  . H FOHCt 

DO  105  J*I.ft  FOHCE 

L*L* 1 FOHCE 

105  P(L)*LOAD(I)*LOAO(J)  FOHCE 

110  CALL  TF.STF  (OFLTA.ACCUK.M.ICNVG)  FOHCt 

120  IF(ICNVG)  400. 130. AGO  FOHCE 

ICNVG  = 0 FAS  NOT  CONVFHGfcl)  FOHCE 

ICNVG*  1 AIL  DELTAS  AH F LESS  Th«N  ACCOx  FOHCE 

130  CALL  GMfcQF (EPSI.tPSO.P)  FOHCE 

IF  ( IT.FO.ITMAX)  I F H=  1 FOHCE 

300  CONTINUE  FOHCE 

400  RETURN  FOHCt 

F.NO  FOHCE 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


SL8RCUT INF  GftPXFA  (A.B.H.ft.N.L) 

PURPOSE 

MIL  TIPL  Y TV  l)  GF  NEPAL  PATH  ICES  TO  Fokm  A 
HESULTANT  GENERAL  MATRIX 

DtSCHIPTICM  OF  PARAMETERS 

A INPLT  FIHS1  MAT  W J X NAftt 
H INPLT  SECCNO  MATHlA  NAMt 
H OUTFtlT  MATRIX  NAftE 

M INPLT  NLMHFh  OF  HCftS  IN  MA(H1A  A OH  H 
N INPLT  NLMPEH  OF  COLUMNS  IN  A Ai.f.  HO«*S  IN  8 
L INPLT  NLMPtR  OF  COLUMNS  IN  MATRIX  b OR  H 

METHOD 

THE  M EY  N MAThIX  A IS  RCSTmuLI IPLIED  BY  ThE  N HY  L MATRIX 
d AND  The  StSULT  IS  STOREU  IN  ThE  M BY  L MATRIX  R . 
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REFERENCE 

NASA  TN  0-6660 


APPENDIX  C (CONT.) 


DIMENSION  A(l)iEll) »Pd> 

lfi»0 

IKs-N 

00  30  K*1*L 

IK*IK*N 

00  20  J*  1 *M 

IR«IR*1 

JI*J-M 

IP*IK 

R(IH)*0. 

DO  10  1*1. N 

JI*JI*M 
I B* I B ♦ 1 

11  R(IR)*R(IR)*A(JI)*H(IB) 
10  CONTINUE 
20  CONTINUE 
30  CONTINUE 
RETURN 
ENU 


SUBROUTINE  GMSHF (A.H.R.MN) 

PURPOSE 

SUBTRACT  ONE  GENFHAL  MATRIX  FRCM  ANOTtit  R TO  FORM 
A RESULTANT  GENERAL  MATRIX 

DESCRIPTION  OF  parameters 

A INPUT  NAME  CF  FIRST  MATRIX 
H INPUT  NAME  CF  SECCNU  MATRIX 
H OUTFUT  MATRIX  NAME 

MN  INPUT  NUMBER  OF  ELEMENTS  IN  MATRIX  A*  H.  OR  R 
METHOD 

EACH  ELFMENT  OF  MATRIX  R IS  SUBTRACTED  FROM  THE 
CORRESPONDING  ELEMENT  OF  MATRIX  A AND  THE  RESULT  IS 
PLACED  IN  THE  CORRESPOND ING  LLEMtNl  UF  MATRIX  R 

REFERENCE 

NASA  TN  0-6660 

DIMENSION  A(1 ) «P ( 1 ) »R<  1 ) 

00  10  IJ=1.MN 
R(IJ)*A(IO)-B(Iw> 

10  CONTINUE 
RFTUHN 
END 


SUBROUTINE  TESTF < A * H « MN « L*  ) 

PURPOSE 

TEST  TRE  APSCLUTF  VALUE  CE  EACE  ELEMENI  CF  MATRIX  A TO 
DETERMINE  IF  IT  IS  LESS  THAN  OR  EQUAL  TO  ThF. 
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APPENDIX  C (CONT.) 

CORRESPONDING  ELEMENT  OF  MATRIX  b 

DESCRIPTION  OF  PARAMETERS 

A INPLT  FIRST  MATRIX  NAME 
B INPLT  SECOND  MATRIX  NAME 

MN  INPLT  NUMBER  OF  ELEMENTS  IN  MATRIX  A CH  b 
LE  OUTPUT  COMPARISON  OF  MATRICES  A ANU  B 

MFTHOD 

IF  ABSOLUTE  VALUE  OF  A(Io)  IS  LESS  THAN  CR  EQUAL  TO 

B ( I J)  FOR  ALL  IJ  = 1*  MN  THEN  LE  « 1 

OTHERWISE.  LE  * 0. 

REFERENCE 

NASA  TN  D-68F0 

DIMENSION  A ( „ ) »E  ( 1 ) 

LF*0 

DO  10  iJslfMN 

IF<ABS(A(IJ) )-B(lJ) ) 10*10, 20 
10  CONTINUE 
LE  *1 

20  RFTURN 
END 


SUBROUTINE  GMEQF < A *R*RN ) 

PURPOSE 

EQUATE  ONE.  GFNFRAL  MATRIX  TO  ANOTHER  GENERAL  MATPIX 

DESCRIPTION  OF  PARAMETERS 
A INPLT  MATRIX  NA«E 
R OUTPUT  MATRIX  NAME 

MN  INPLT  NUMBER  OF  ELEMENTS  IN  MATRIX  A CR  H 
METHOD 

EACH  ELEMENT  OF  MATRIX  R IS  SET  FOUaL  TO  ThE 
CORRESPONDING  ELEMFNT  OF  MATRIX  A 

REFERENCE 

NASA  TN  0-6860 

DIMENSION  A ( 1 ) « R ( 1 ) 

DO  10  I J* 1 » MN 
R(IJ)xAIL) 

10  CONTINUE 
RFTURN 
END 
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